The flow induced by a mock whale : origin of flukeprint formation. 
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Fluke oscillations of cetaceans (whales, dolphins and porpoises) imprint the ocean surface with 
an oval patch that surface waves cannot enter. Scientists have hypothesized that modification 
of surface tension by surfactants creates flukeprints. Here, we show, on the contrary, that the 
formation of flukeprint is primarily due to a wave-current interaction problem. We provide Particle 
Image Velocimetry measurements of anisotropy and vorticity in the flow generated by a mock whale 
fluke in a laboratory experiment. We explain for the first time why long gravity waves enter the 
flukeprint, whereas short gravity waves are blocked. 



PACS numbers: 47.63.-b, 47.63.M-, 47.63.mc 

Free propagation of waves is an exception. For exam- 
ple, in a marine environment, water waves are refracted 
by the varying landscape of the ocean floor as they ap- 
proach the seashore and ocean currents modify wave 
propagation. We consider wave interactions that have 
fascinated both casual whale-watchers and marine biol- 
ogists: cetacean flukeprints. Cetacean flukeprints (made 
by whales, dolphins and porpoises) occur when a vortex 
ring, shed by the fluke, interacts with waves present on 
the free surface of the water [1 . In the ocean these pat- 
terns are sometimes known as "whale footprints". The 
surface signatures are characterized by a smooth oval pat- 
tern where short gravity waves cannot penetrate. De- 
pending on flow conditions, one can also observe several 
radii surrounding the smooth oval corresponding to mode 
conversion into blue-shifted and capillary waves as can be 
seen for dolphins flukeprints easily observed in a swim- 
ming pool (Fig. [T]). The key explanation of flukeprint 
formation, first revealed in the present work, involves the 
dispersive properties of water waves along with the effect 
of gravity, as well as surface tension and wave-breaking 
in the mode conversion induced by the surface current. 

Similar surface patterns occur when a vortex ring cre- 
ated by an upward-pointing jet interacts with a free sur- 
face, although these patterns are simpler due to the reg- 
ular shape of the jet orifice. Another related example 
is a circular hydraulic jump in which a jet of water im- 
pacts a solid plate and creates a disk-shaped region with a 
boundary that surface waves cannot penetrate [2 . Long 
shallow water waves are blocked by the flow induced on 
the plate by the impacting jet with a characteristic bor- 
der, namely the jump. In essence, the flukeprint is an 
inverted circular jump caused by a jet-like flow induced 
by a vortex ring obliquely impacting the ocean surface 
with different dispersive properties. 

Here, we provide for the first time, PIV measurements 
from experiments reported in [1 where an oscillating 2D 
polypropylene fluke created vortex rings and resulting 
surface prints were visualized by particles and a laser 
sheet. 

The landscape of waves on the ocean surface is com- 
plex, with interactions of both gravity and capillary 
waves [7 . As argued in [1 , the current forming a 




FIG. 1: Dolphin's flukeprints. Top: A close-up of the edge of a 
dolphin flukeprint where mode-converted capillary waves can 
be observed near the boundary and are subsequently damp- 
ened in the print. Bottom: A sequence of four flukeprints 
created by a dolphin. The flukeprints grow in size over time, 
so the small prints are the ones made most recently by the dol- 
phin traveling from "bottom to top" (dolphin fluke is visible 
at the top of the picture). 



flukeprint originates below the surface at the swimming 
depth of the whale. As a vortex ring forms, the current 
flows upwards through the center and then radially out- 
wards at the ocean's surface, interacting with oncoming 
gravity and capillary waves already present. 
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In deep water and without surface tension, when 
a gravity wave meets a counter-current, the incident 
wavelength diminishes (the phenomenon known as blue- 
shifting) because of normal dispersion. That is, the wave 
group velocity is a decreasing function of the wavenum- 
ber k. In this regime {kh » 1 where h is the water 
depth), the blocking of pure gravity waves (neglecting 
surface tension 7 = 0) is then described by the following 
formula [4|: 



where T is the period of the incoming waves in the gravity 
field g. This equation indicates that the blocking velocity 
is linearly related to the wave period. We expect wave 
blocking if the period of the existing gravity waves is 
short, but no blocking if the period is long. This simple 
formula explains the presence of only long waves inside 
the flukeprints and has not previously been noted. Since 
the blocking velocity in deep water depends on the wave 
period, it is an intrinsically dispersive effect whereas in 
shallow water the blocking occurs for a velocity which 
depends only on the water depth h that is Uh = \/~gh 

mm- 

When including surface tension (7 ^ 0) in the grav- 
ity dispersion relation, Badulin et al. [11 observed that 
as long as the counter-current is strong enough, grav- 
ity waves can still be blocked in deep waters. In addi- 
tion, blue-shifted waves appear by mode conversion but 
are also stopped while drifting backward at a new block- 
ing boundary, formed where the blue-shifted wave merges 
with a new capillary solution O [TTHT3] . Capillary waves 
appear at this secondary blocking line, described in [5], 
and propagate in the same direction as the incident grav- 
ity waves, provided the regime stays linear. Ultimately, 
the capillary waves no longer propagate far inside the 
flukeprint because they are not only quickly damped by 
viscosity but they are no longer solutions of the dispersion 
relation [5 , and are seen on its boundary. The threshold 
velocity for the appearance of capillary waves in deep wa- 

/ xl/4 

ter on a flow current is [/^ — 27r ^ where = >/2 ( ] 

is the minimum of the phase velocity as a function of the 

wavenumber and where Ic = 'sj^ is the capillary length 

[5]. If the amplitude grows too much, wave breaking 
will lessen the conversion of the incident gravity waves 
to capillary waves. In this case, energy is dissipated in 
turbulence rather than transferred linearly towards the 
capillary waves. 

In addition to its influence on the mode conversion, 
surface tension can play another role in the flukeprint for- 
mation since it also changes because of variation in sur- 
factant concentration or temperature. This phenomenon 
in which surface tension is non-increasing as the temper- 
ature or surfactant concentration increases is referred to 
generically as a Marangoni effect. Oily skin sloughed by 
cetaceans may also modify the properties of the air- water 



interface, but in a different way than just changing the 
value of surface tension. An oil monolayer can feature a 
surfactant concentration gradient (Gibbs-Marangoni ef- 
fect) because of the motion of water waves [14 . This 
induces a surface dilatational elasticity (Gibbs surface 
elasticity) related to the presence of elastic longitudinal 
waves within the monolayer. The resulting rigidity of 
the surface dampens the capillary waves (which consti- 
tutes the Franklin effect) by increasing the shear forces 
on the free surface hence the viscous dissipation. More- 
over, capillary waves are naturally dampened by viscosity 
and their amplitude is known to decreases exponentially 
with a factor proportional to the water viscosity /i and in- 
versely proportional to surface tension 7 and wave period 

27r/x 

T: A = Aqc [14]. Hence, a decrease of surface ten- 
sion quickly dampens capillary waves. In addition to the 
blocking of short gravity waves by the induced flow cre- 
ated by the vortex ring, the Franklin effect may provide 
another mechanism for the damping of capillary waves in 
addition to their usual dissipation by viscosity. 

Recently, it was demonstrated [8 using thermal im- 
agery that the inside of flukeprint has a colder water tem- 
perature than the outside. This temperature signature is 
a natural consequence of the vortex ring inducing a jet of 
water taking cold water from below the free surface and 
dispatching it at the surface forming the interior of the 
flukeprint. The thermal signature provides one way to 
track whales. However, it must be noted that the ther- 
mal change is not the origin of the print. A decrease 
in temperature will increase the surface tension, hence 
decrease the viscous damping of the capillary waves. 

To visualize vortex ring formation and the resulting 
surface patterns. Levy et al. P conducted a series of ex- 
periments in collaboration with Flo-Metrics in San Diego 
using a tank of water approximately 122 cm long by 61 
cm wide by 49 cm deep. In water, they oscillated a 2D 
fluke shape made of slightly flexible polypropylene (0.16 
cm thick) with the boundary shape scaled up from a pho- 
tograph of a blue whale fluke. The mock fluke was 0.2 
m tip to tip, whereas blue whale flukes are typically 5-6 
m. The velocity field was computed by particle image 
velocimetry (PIV) using a laser sheet and small parti- 
cles similar to Kaleiroscope. The laser sheet illuminating 
the small flakes was created using a 1 Watt green laser 
directed through a 3 mm diameter stirring rod. 

In each experimental run, oscillations were maintained 
at 0.32 Hz and no forward motion. For comparison, blue 
whale fluke oscillations occur at about 0.5-1.0 Hz with 
forward velocity 1-1.5 m/s. The initial angle of the fluke 
with the plane of the surface was varied, and recorded 
video at 30 fps. Stills from this video were first reported 
in [1 but did not include PIV. Here, PIV measurements 
from several runs initiated at 12.5 and 60 degrees from 
horizontal are reported. Each trial consisted of one os- 
cillation from the prescribed angle, down 25 degrees and 
then up. While the motion of the robotic fluke in the 
experiments is only a rough approximation of whale mo- 
tion, the observations are consistent with the theory de- 
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FIG. 2: Surface flows. Top: Typical instantaneous velocity 
field on the free surface created by an impacting vortex ring 
(t=0.266s after flapping begins). The fluke is oriented as if 
the whale is facing left and the flow is stronger on average on 
the right. The colors correspond to the modulus of the veloc- 
ity at each point. The angle of the fluke at the maximum is 
60° and it oscillates through a 25° arc. Bottom: Correspond- 
ing strikelines (perpendicular to the streamlines of the flow) 
which display an oval-shape flow print similar to the surface 
flukeprint as seen by the observers. The colors correspond to 
the modulus of the velocity at each point. 



water from the center of a patch with an oval shape, see 
Fig. [2] The induced surface current is anisotropic and 
time-dependent; it is stronger on the rear of the flukeprint 
than on the front. Our Fig. [2] (top) is very similar to the 
numerical simulations reported in 17 where a vortex 
ring impacts obliquely a free surface. On average, the 
flow is weaker inside the flukeprint than on its boundary. 
The streamlines (not reported here) features a null-line 
from which all streamlines radiate outward. The strike- 
lines (perpendicular to the streamlines) display a typical 
oval shape reminiscent of field observations, see Fig. [2] 
(bottom). 




scribed in the previous section. The image processing 
is completed on a personal computer with the software 
Davis 7.2 from Lavision. Coordinate x denotes the cruise 
direction of the whale, y the transverse direction and z 
the vertical direction. 

According to Archer et al. [15 , the life-time of a vor- 
tex ring interacting with a free surface at 90° of incidence 
features three stages after creation: approach of the free 
surface, slowing and expansion. The PIV measurements 
show that an outgoing flow is induced on the surface of 



FIG. 3: Vortex rings. Top: Back view of the whale (t=4.2s 
after flapping begins). Typical quadrupolar flow on the rear 
of the fluke. Cross- sections of the two vortex rings. The angle 
of the fluke is 60°. Bottom: Side view (t=3.6s after flapping 
begins). Streamlines of the flow. One observes the cross- 
section of the two vortex rings created during the downward 
and the upward motions of the fluke. The upper ring creates 
free surface shear flows which will block the waves. The colors 
correspond to the modulus of the velocity at each point. The 
angle of the fluke is 60°. 
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In Fig. [3] (top) a quadrupole is seen behind the mock 
fluke. The quadrupole corresponds to the successive 
emissions of two vortex rings observed from the side: one 
towards the deep in the downward motion of the fluke, 
one towards the free surface in the upward motion of the 
fluke, see Fig. [s] (bottom). Measurements of the surface 
flow profile in the radial direction demonstrate the in- 
homogeneity of the flow as one gets far from the center 
(y = 0) of the flukeprint, see Fig. |4] The radial veloc- 
ity increases almost linearly with the distance to center 
before reaching a maximum and then decreases to zero. 

The radial flow depends on the vertical position as seen 
on the Fig. [4] (bottom). In these experiments, the PIV 
indicates that the flow is almost linear with the depth 
with a flow reversal below the core of the toroidal vortex 
tube. 

As a conclusion, the existing theory for wave-current 
interactions includes a number of simplifying assump- 
tions. The first, of course, is that we assume the flow 
is in one or two dimensions rather than providing a fully 
3D model. Even in the 2D case, our measurements show 
that the surface flow generated in the whale flukeprint 
is anisotropic in the radial direction. Another simplifica- 
tion is the assumption that the flow has usually a uniform 
velocity profile on the entire depth, or on a given depth 
[6^ or at most a linear velocity profile on the entire depth 
[l8.-.26j or on given depth [6 . The PIV displays a flow re- 
versal in the vertical direction in addition to a maximum 
in the radial direction. 
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FIG. 4: Velocity profiles. Top: Back view of the whale (z=- 
10mm below the free surface). Typical longitudinal variation 
of the radial surface flow induced by the upper vortex ring 
(t=2.6s after flapping). Bottom: Back view of the whale (x=- 
132mm on the left of the upper vortex ring center). Typical 
vertical variation of the radial surface flow induced by the up- 
per vortex ring (t=8.4s after flapping). Starting at an initial 
12.5° from the surface, the fluke oscillates through an arc of 
±25°. 



